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Solving for closure errors due to polarization leakage in radio 
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Abstract. Mechanical and electronic imperfections can result into polarization leakage in individual antennas of a 
radio interferometer. Such leakages manifest themselves as closure errors even in co-polar visibility measurements 
of unpolarized sources. This paper describes and tests a method for the study of polarization leakage for radio 
interferometric telescopes using only the nominally co-polar visibilities for unpolarized calibrators. Interpretation 
of the resulting closure phases on the Poincare sphere is presented. Since unpolarized sources are used, the actual 
solutions for leakage parameters is subject to a degeneracy which is discussed. This however, does not affect the 
correction of closure errors in our scheme. 
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1. Introduction 



Stokes coherence function for a source of known structure 



The m jumal coherence ftrrottan (also called Lhc vlslbil- ct aL ^ henceforth HBS) 



(usually an unresolved source) ( Hamaker et al. 1996 ; Bault 



ity function) for an unresolved and nnprferfcert source, 

measured by an interferometer array can be modeled as a 
product of antenna based complex gains. These complex 
gains can be derived from the measured visibility func- 
tion using the standard algorithm, which we call antsol. 
antsol forms the central engine of most amplitude and 
phase calibration schemes used for radio interferometric 
data. (The earliest published reference for an algorithm for 
antsol of which we are aware is Thompson & D'Addario 

Usually one measures the components of the incident 
radiation along two orthogonal polarization states by us- 
ing two separate feeds. The signals from the two feeds 
travel through essentially independent paths till the cor- 
relator. However, due to mechanical imperfections in the 
feed or imperfections in the electronics, the two signals 
can leak into each other at various points in the signal 
chain. 

At the correlator, signals from all the antennas are 
multiplied with each other and the results averaged to 
produce the visibilities. The signals of same polarization 
are multiplied to produce the co-polar visibilities while the 
signals of orthogonal polarizations are multiplied to pro- 
duce the cross-polar visibilities. The co-polar and cross- 
polar visibilities can be used to compute the full Stokes 
visibility function. Antenna based instrumental polariza- 
tion and polarization leakage can be derived from the full 
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The correlator used for the Giant Metrewave Radio 
Telescope (GMRT) by default computes the co-polar visi- 
bilities using the Indian mode of the VLB A Multiplier and 
Accumulator (MAC) chip. Here we describe a method, 
which we call leaky antsol, for the computation of the 
leakages using only the co-polar visibility function for an 
unpolarized source. Following the notation used by HBS, 
we label the two orthogonal polarizations by p and q to re- 
mind us that the formulation is independent of the precise 
orthogonal pair of polarization states chosen. 

Sect. H describes the motivation which led to this anal- 
ysis. For orientation, Sect. || starts with the problem of 
solving for the usual complex antenna based gains and 
sets up an iterative scheme for the solution. The problem 
of simultaneously solving for the complex antenna gains 
and leakages is then posed in Sect. 3.1 and a similar it- 
erative scheme is set up. Sect. 3J2 presents the results of 
the simulations done to test the scheme. Sect, hi presents 
some results using the GMRT at 150 MHz. Also, we were 
fortunate to have the L-band feeds of one of the GMRT an- 
tennas converted from linear to circular polarization. We 
observed 3C147 in this mode where all baselines with this 
special antenna measured the correlation between nomi- 
nally linear and circular polarization. Results of this ex- 
periment demonstrate that the leakage solutions are in- 
deed giving information about the polarization proper- 
ties of the feeds. These results and their inter pret ation on 
the Poincare sphere are presented in section h2. Sect. || 
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gives the interpretation of the leakage solutions and dis- 
cusses closure errors due to polarization leakage using the 
Poincarc sphere. 

2. Motivation 

Rogers (1983) pointed out in the context of the VLBA, 
that the polarization leakage cause closure errors even in 



nominally co-polar visibilities. Massi et al. (1997) have car- 
ried out a detailed study of this effect for the telescopes of 
the European VLBI Network (EVN). Our motivation in 
this paper is that the current single sideband GMRT cor- 
relator uses the so called Indian mode of the VLBA MAC 
chips to produce only the co-polar visibilities. Also, the 
planned Walsh switching has not yet been implemented 
at the GMRT and in any case, would not eliminate leak- 
age generated before the switching point. Tests done using 
strong point source dominated fields show unaccounted 
closure errors at a few percent level. The motivation be- 
hind developing an algorithm to solve for gains and leak- 
ages simultaneously, using only the co-polar visibilities was 
to determine if the measured closure errors could be due 
to polarization leakage in the system. Estimates of leak- 
age can then be used in the primary calibration to remove 
the effects of polarization leakage. This is where this paper 
differs from the earlier work of HBS which is about the cal- 
ibration using the full Stokes visibility function, needed for 
observations of polarized sources. The polarization leak- 
age in some of the EVN antennas corrupts the co-polar 
visibilities at a level visible as a reduction in the dynamic 
range of the maps (Massi & Aaron 1997a ^; Massi ct al 



199S). Thus such a method can also be used in imaging 



data from the EVN and other telescopes affected by such 
closure errors. 

Let gf represent the complex gain for the p- 
polarization channel of the i th antenna and a q repre- 
sent the leakage of the (^polarization signal into the p- 
polarization channel. The electric field measured by an- 
tenna i can then be written as 



E p = q p E p +a q E q 

t Ji 1.0 I "» ^l.O 



(1) 



where Ef and E? Q are the responses of an ideal an- 
tenna to the incident radiation in the p- and g-polarization 
states respectively^]. For an unpolarized source of radi- 
ation, (Ef Ej* ) = 0. The co-polar visibility for such a 
source, measured by an interferometer using two anten- 
nas denoted by the subscripts i and j, is given by 



9z9i 9 



q q* qq 
*-*7 Pij.c 



(2) 



1 No te that Eq. [I] is equivalent to the Eq. 1 of Massi et al 
(1997 ), who use a different parameterization. Strictly speaking 
the ratio af/gf is the correct measure of the leakage and this 
is what has been plotted in FigJ^ below. This does not however 
affect the computational scheme described here. 



where is independent gaussian random baseline 
based noise and p pp Q = (Ef E p * ) and pff a = (Ef Ej*) 
are the two ideal co-polar visibilities, usually represents 
the contribution to p pp which cannot be separated into an- 
tenna based quantities, therefore is a measure of the 
intrinsic closure errors in the system and is usually small. 

For an unpolarized point source (Ef E p * a ) = 
(E? a E!j*) = p pp a = 1/2 where I is the total flux density. 
Writing Xf p = pfj / p PP Q we get 



X PP 



p p* 



i Q Q* 



+ e 



ij 



(3) 



where tij now refers to the baseline based noise in X pp . 
Assuming a's to be negligible, the usual antsol al- 



gorithm estimates gfs such that ^2 



VPP 

\ x iJ 



V p + 



is 



ed (see section ||). Normally, Walsh switching 



(Fhompson et al. 1986) is used to eliminate the polar- 
ization leakage due to cross-talk between the signal paths, 
such that aja q * <C eij. However, afs can also be finite 
due to mechanical imperfections in the feed or the cross- 
polar primary beam, which cannot be eliminated by Walsh 
switching. 

In the case of significant antenna based polarization 
leakage (compared to JTij) , the second term in Eq. % in- 



volving a\ fs will combine with the closure noise e»j. The 
polarization leakage therefore manifests itself as increased 
closure errors (see Sect. || for a geometric explanation on 
the Poincare sphere). 



3. Algorithm and simulation 

In the absence of any polarization leakage, g^s can be es- 
timated by minimizing 



El* 



-9i9 P K 



(4) 



with respect to <jfjS, where w pp = l/<7jy, 0y being the 



variance on the measurement of X 
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In Eq. |2|, if p PP accurately represents the source struc- 



ture, Xf p will have no source structure dependent terms 
and is purely a product of two antenna dependent com- 
plex gains. For a resolved source, p PP Q can be estimated 
from the image of the source. 

Evaluating t^It and equating it to zerop|( pnatnagar 



199J), we get 



T X pp q p w pp 



p j^i 



(5) 



This can also be derived by equating the partial deriva- 
tives of S with respect to real and imaginary parts of gf . 

2 Complex derivatives can be evaluated by treating gf and 
gf* as independent variables (Palka 199C). 



Bhatnagar & Nityananda: Measurement of polarization leakage 



3 



Since the antenna dependent complex gains also ap- 
pear on the right-hand side of Eq. H, it has to be solved 
iteratively starting with some initial guess for gjs or ini- 
tializing them all to 1. Eq. can be written in the iterative 
form as: 



91 + x 



E 

j 



(6) 



where n is the iteration number and < A < 1. 
Convergence would be defined by the constraint 
\S n — S n -i | < (3 (the change in S from one iteration to an- 
other) where, (3 is the tolerance limit and must be related 
to the average value of tij . Eq. || forms the central engine 
for the classical antsol algorithm used for primary cali- 
bration of the visibilities and in self-calibration for imaging 



purposes. This algorithm was suggested by Thompson & 



D'Addirio (1982) 



3.1. The leaky antsol 

In the presence of significant polarization leakage, Eq. [| 
can be used to re-write Eq. 



as 



5': 



E 



yPP 



(gig 



a ! a i ) 



w- 



(7) 



In this form 
noise e 



S is an estimator for the true closure 



,j rather than the artificially increased closure noise 
a i a j e *J') °- uc ^° ^ nc P resence °f polarization leakage. 
Equating the partial derivatives to zero, we 



get 



<l 



j 

p s¥i 



p w pp 
3 y 



ofl a 9 q p w PP 



E \gS\ 



(8) 



w,, 



self-calibration, when only relative phases of antennas are 
determinate, with one antenna acting as an arbitrary ref- 
erence. For observations of unpolarized sources, we can 
similarly say that any feed can be chosen as a reference 
polarization, with zero leakage, and other feeds have gains 
and leakages in the basis defined by this reference. Other 
conventions may be more convenient, as discussed in the 
appendix which discusses degeneracy in detail. 

3.2. Results of the simulations 

We simulated visibilities with varying fraction of polar- 
ization leakage in the antennas to test the algorithm as 
follows. The antenna based signal and leakage were con- 
structed as gi = R g and on = f ■ R a where R g and 
R a were drawn from the same gaussian random popula- 
tion. The visibility from two antennas i and j was then 
constructed as Xij = gig* + a^a^ + Cij for < / < 0.1. 
This is equivalent to a visibility of an unpolarized point 
source of unit strength with a complex antenna based gain 
gi and leakage ai of strength proportional to /. Eq. [6] was 
then used to compute gi 

2 



and residual x 2 computed as 



i 



SiS* 



The computed values of gf were 

then used to compute improved estimates for g p by simul- 
taneously solving for gf and a? using the iterative forms 
of Eq. H and |j. The derived values of gf and af matched 
the true values to within the tolerance limit. A new % 2 

2 



1 



The val- 



was computed as yf = V\- „• . 

ues of Xa an d X\ as a function of / are plotted in Fig. |[ 
The two curves become distinguishable when the leakage 
is significantly greater than (for / greater than ~ 1%). 
After that, the value of X\ is consistently lower than Xai 
where the contribution of antenna based leakage has not 
been removed. Also notice that X\ remains constant while 
X 2 quadratically increases as a function of /. This is due 
to the fact that antsol treats the antenna based polariza- 
tion leakage as closure errors resulting in an increased x 2 
with increasing fractional leakage. 



4. Real data 



E Xffalwff - gf E 9i*a q iW: 



VP 
ij 



E 



(9) 



These non-linear equations can also be iteratively 
solved. 

Eq. |^, which expresses the observed visibilities on a 
point source unpolarized calibrator in terms of the gains 
and leakage coefficients of the antennas, would take the 
same form if written in an arbitrary orthogonal basis. It is 
clear that the g's and the a's will change when we change 
the basis, so this means that the equations cannot have a 
unique solution. This situation is familiar from ordinary 



4.1. 150 MHz data 

Engineering measurements for polarization isolation at 
150 MHz for the GMRT show significant polarization leak- 
age in the system. We therefore used leaky antsol to cal- 
ibrate the data from the Galactic plane phase calibrator 
1830-36 which is known to be less than 0.2% polarized at 
1.4 GHz. The percentage polarization at 150 MHz is not 
known, but it is expected to decrease further and it was 
taken to be an unpolarized point source. 

Fractional polarization leakage (\af/gf |) of up to 100% 
was measured for most of the antennas, which is consis- 
tent with the estimated leakage measured from system 
engineering tests. Again, XaL ana - X\ were computed and 
the results are shown in Fig. |. The 150-MHz GMRT band 
suffers from severe radio frequency interference (RFI). The 
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Fig. 1. Figure showing the results of the simulations. The 
top curve is the value of x 2 using the classical antsol 
(Xa)- The bottom curve is the value of x 2 using the 
leaky antsol {x\) as a function of the percentage po- 
larization leakage. 
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Fig. 2. Figure showing the results using GMRT data at 
150 MHz for the compact Galactic plane source 1830-36. 
The top curve is the value of x 2 using the classical antsol 
(Xa)- The bottom curve is the value of x 2 using the 
leaky antsol {x 2 ) as a function of time. 



sharp rise in the value of x\ around sample number 10 is 
due to one such RFI spike. This spike is present in the 
total power data from all antennas at this time. On an 
average, the x 2 reduces by ~ 60% when leakage calibra- 
tion is applied (xi)- This is consistent with polarization 
leakage being a major source of non-closure at this fre- 
quency. 

4.2. L-band data with circular feed 

The GMRT L-band feeds are linearly polarized. For the 
purpose of a VLBI experiment conducted in December 
2000, the L-band feed of one of the antennas was con- 
verted to a circularly polarized feed. The rest of the L- 
band feeds were linearly polarized and we took this op- 
portunity to measure correlations between the circularly 
polarized antenna with other linearly polarized antennas 
using the source 3C147. Two scans of approximately one 
hour long observations were done using the single side 
band GMRT correlator. This correlator computes only co- 
polar visibilities. With this configuration of feeds, visibili- 
ties between the circularly polarized antenna and all other 
linearly polarized antennas corresponds to correlation be- 
tween the nominal X- and R-polarizations, labeled by RX, 
were recorded in the first scan. The polarization of the cir- 
cularly polarized antenna was then flipped for the second 
scan to record the correlation between the nominal X- and 
L-polarization states, labeled by LX. 



The VLA Calibrator Manual]^] lists the percentage po- 
larization ( ^ Q2+ f 2+V2 ) for 3C147 at L-band < 0.1%. The 
cross-polar terms in Eq. |^, which are assumed to be zero, 
will therefore contribute an error of the order of 0.1%. 
These cross-polar terms will be, however, multiplied by 
gains of type gf af . Since gf and a\ are both assumed 
to be uncorrelated between antennas, this error will man- 
ifest as random noise in Eq. ||. Within the limits of other 
sources of errors, the source 3C147 can therefore be con- 
sidered to be a completely unpolarized source. 

4.2.1. Results and Interpretation 

The response of an ideal circularly polarized antenna to 
unpolarized incident radiation can be expressed as a su- 
perposition of two linear polarization states as Ef = 

Ei C c , e '' S + EY.o e ~ iS where, the superscripts R, X and Y 
denote the right circular and the two linear polarization 
states respectively. S is half the phase difference between 
the two linear polarization states and is equal to 7r/4 for 
right-circular polarization and — 7r/4 for left-circular po- 
larization. Writing the general Eq. ^ for right-circularly 
polarized antenna as Ef- = gfEf Q + a\E^ and substi- 
tuting for Ef Q and E^ a we get 



The VLA Calibrator manual is available on the web from 



tittp : //www . aoc .nrao . edu/- gtaylor/ calib .html 
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Fig. 3. Figure showing the results using visibilities with 
one circularly polarized antenna and all other linearly po- 
larized antennas at L-band. The x- and y-axis denote the 
real and imaginary parts of af / gf respectively, gf and a\ 
were solved for every integration time (~ 17s). All linearly 
polarized antennas are close to the origin in this plot. The 
solutions for the circularly polarized antenna (C03) are 
the set of points away from the origin (shown by open 
circles and triangles) . The two sets of points for this an- 
tenna, separated from each other by ~ 180° are solutions 
for the right- and left-circular polarization channel. The 
points denoted by open circles are from correlation be- 
tween the right- and left-circular polarization of C03 with 
nominal linear X-polarization of the other antennas (la- 
beled as C03 RX and C03 LX respectively). The points de- 
noted by triangles are from correlation of C03 with nomi- 
nal Y-polarization of the other antennas (labeled as C03 RY 
and C03 LY ). 



E? =gf (E^ S + Ej^) + 

of (EX e-" + E% e*) ' ' 

Eq. H for the case of correlation between a circularly 
polarized and a linearly polarized antenna, with polariza- 
tion leakage in both the antennas, can be written as 

X% x =(g^ s + a^ s )gf + 

(gfe~' s + afe^)af (11) 

_ X' X* . Y' Y* 

=9i 9j + a i <Xj 

where gf = gfe LS +afe~ LS and af = gf e~ a +af e lS . 
The leaky antsol solutions for the circularly polarized 
antenna in this case will correspond to gf and aj . 



Let Pi — of I gf (Pj = af j gf for the circularly po- 
larized antenna). Then, the amplitude of Pi is a measure 
of the fractional polarization leakage in the antenna while 
the phase of p gives the phase difference between the 
signal from one of the feeds and the leaked signal from 
the other feed. For an ideal circularly polarized antenna, 
| Pi | w 1. A plot of the real and imaginary parts of this 
quantity for all antennas should therefore clearly show Pj 
for the circularly polarized antenna with an amplitude of 1 
and at an angle of tt/2 with respect to the nominal X-axis. 

The real and imaginary parts of Pi for all antennas 
from this experiment are shown in Fig. ||. The solutions 
were computed for every integration cycle of ~ 17 sec 
and the points on this plot represent the tip of phasor 
Pj . The collection of points near the origin are for all the 
linearly polarized antennas while the collection of two sets 
of points away from the origin, approximately an angle of 
7r from each other, are for the circularly polarized antenna. 
The solutions found by leaky antsol match the expected 
results quite well. This therefore constitutes a reasonably 
controlled test with real data showing that the solutions 
indeed provide information about the polarization leakage 
in the system. 

This experiment however provides much more infor- 
mation about the polarization properties of the antenna 
feeds used. The collection of points in the first quadrant 
denoted by open circles are the values of Pi derived from 
the correlation between the nominal right-circularly po- 
larized signal and the linearly polarized signals along the 
nominal X-axis from all other antennas. Points in the third 
quadrant are similarly derived using the left-circular sig- 
nals. The set of points denoted by triangles in the sec- 
ond and fourth quadrant are derived using correlations of 
right- and left-circularly polarized signals with the linearly 
polarized signals along the nominal Y-axis from all other 
antennas. 

A larger spread in the solutions using the left-circularly 
polarized signals indicates that the closure noise (from 
other unknown sources) in these signals is higher. The fact 
that the amplitude of Pi derived using the right-circularly 
polarized signals is rs 0.5 indicates that the nominal cir- 
cularly polarized feed is in fact elliptically polarized with 
this axial ratio. The spread of ±1 — 2% about the ori- 
gin is indicative of polarization leakage at the level of few 
percent in the linearly polarized antennas as well. The 
leakage in one of the linearly polarized antennas is signif- 
icantly larger (w 4%). Since this kind of data is routinely 
taken on primary calibrators during GMRT observations 
for synthesis imaging, leaky antsol provides a useful di- 
agnostic of system health, polarization performance and 
numbers needed to correct the data in high accuracy work. 

The following test was also carried out to check that 
the closure phase on a triangle involving the circular feed 
was indeed mainly due to polarization effects. The three 
baselines making up this triangle were flagged as bad base- 
lines from the input data and a new solution found for the 
gains and leakages of all antennas. This solution predicted 
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resented by the vector 



i<t> 



Clearly, 9 = n/2 



the same closure phase (to within errors) as actually ob- 
served. 



5. Closure phase and the Poincare sphere 

In this section we use right- and left-circular polariza- 
tion states as the basis. A general elliptically polarized 
state can be written as a superposition of two states rep- 

cos6>/2 
sin6>/2 e- 

corresponds to linear polarization and 9 ^ 0, 7r/2 to ellip- 
tical polarization. Increasing <fi by £ rotates the direction 
of the linear state or the major axis of the ellipse by £/2. 
We can chose the phase of the basis so that = cor- 
responding to linear polarization along the x-axis. The 
Poincare sphere representation of the state of polariza- 
tion maps the general elliptic state to the point (9, <p) on 
the sphere. The properties of this representation are re- 
viewed by Ramachandran fc Ramaseshan (1961 ). We are 
concerned here with one remarkable property, discovered 
by Pancharatnam (1956| , 1975| ). Whenever there is con- 
structive interference between two sources of radiation, 
it is natural to regard them as in phase. An unexpected 
property of this simple definition manifests itself when we 
consider 3 sources of radiation of different polarization - 
that if a source A is in phase with B and B in phase 
with C, C in general need not be in phase with A. The 
phase difference between A and C is known in the optics 
literature as the geometric or Pancharatanam phase (see 
also Ramaseshan & Nityananda (1986); Berry (1987)). We 
show that this naturally occurs in radio interferometry of 
an unpolarized source with three antennas of different po- 
larizations. 

Let the polarization states of the three antennas be 



represented by 
lar basis. Denoting the vector 



and 



93 

013 



m a circu- 



by ipi, one clearly 



9i \ I 92 

Oil J ' V «2 , 

' 9i 

Q-i 

see that the visibility on the 1-2 baseline is propor- 
tional to i>iip2- Hence the closure phase around a tri- 
angle made by antennas 1, 2, and 3 is the phase of the 
complex number (also called the triple product) VI23 = 
(V4V , 2)(V4V'3)('03' ! /'i)- In the quantum mechanical litera- 
ture, this type of quantity goes by the name of Bargmann's 
invariant and its connection to the geometric phase was 
made clear by Samuel & Bhandari (1988). With some 



work, one can give a general proof that the closure phase 
(phase of V123) is equal to half the solid angle subtended 
at the centre of the Poincare sphere by the points repre- 
sented by ip\, ip2, and ip3 on the surface of the sphere. For 
the case where the polarization state of the three anten- 
nas are same, this phase is zero in general. However, when 
the polarization states of the antennas are different, this 
phase is non-zero. 

The well known result that an arbitrary polarization 
state can be represented as a superposition of two orthogo- 
nal polarization states translates to representing any point 
on the Poincare sphere by the superposition of two dia- 



metrically opposite states on a great circle passing through 
that point. For example, circular polarization can be ex- 
pressed by two linear polarizations, each with intensity 
1 / v2- In the context of the present work, the nominally 
circularly polarized antenna maps to a point away from 
the equator on the Poincare sphere (it would be exactly 
on the pole if it is purely circular) while the rest of the 
antennas map close to the equator (they would be exactly 
on the equator if they are purely linear and map to a sin- 
gle point if they were also identical). The visibility phase 
due to the extra baseline based term in Eq. |^ due to polar- 
ization mis-match is a consequence of the Pancharatanam 
phase mentioned above. This phase, on a triangle involv- 
ing the circularly polarized antenna, will be close to the 
angle between the two linear antennas. For example, if 
1 \ / 1 \ , , I cos 7 



ipi = [l),^= I q J > and V>3 = sin7 J > tnc Ph ase 

of V123 will be 7. This picture can be depicted by plotting 
the real and imaginary parts of af/gf, which is done in 
Fig. |3|. The circularly polarized antenna can be clearly lo- 
cated in this figure as the set of point away from the origin 
while the linearly polarized antennas as the set of points 
close to the origin. The collection of points located away 
but almost symmetrically about the origin represents the 
nominal right- and left-circularly polarized feeds. Points 
on the equator, but significantly away from the origin rep- 
resents an imperfect linearly polarized antenna. Note that 
the average closure phase between the nominally linear an- 
tennas is close to zero, which defines the mean reference 
frame in Fig. [|. 

6. Conclusions 

Rogers (1983) pointed out that non ideal feed polariza- 
tions of the individual antennas of a radio interferometer 
can result into closure errors in the co-polar visibilities. In 
this paper we describe and demonstrate a method to mea- 
sure the polarization leakage of individual antennas using 
the nominally co-polar visibilities for an unpolarized cali- 
brator. This method can therefore be used as a useful tool 
for studying the polarization purity of the antennas of ra- 
dio interferometers from the observations of unpolarized 
calibrators. However, since only unpolarized calibrators 
are used, the actual solution for the leakage parameters is 
subject to a degeneracy. This degeneracy does not affect 
the correction of the visibilities and can be used to remove 



the closure errors due to polarization leakage. Massi et al _ 
|(1997| ) have shown that such polarization leakage induced 
closure errors in the data from the EVN is the dominant 
effect of instrumental polarization. For the EVN, this ef- 
fect can be seen as a reduction in the dynamic range of the 
images. Our method can be used for such data to remove 
these closure errors for unpolarized sources. 

The general elliptic state of the polarization of radia- 
tion can be represented by a point on the Poincare sphere. 
The phase difference between three coherent sources of ra- 
diation but with different states of polarization goes by the 
name of Pancharatanam or geometric phase in the optics 
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literature. We interpret the co-polar visibilities with po- 
larization leakages on the Poincare sphere and show that 
the polarization induced closure phase errors in radio in- 
terferometers is same as the Pancharatanam phase of op- 
tics. The antenna based leakages also map to points on 
the Poincare sphere and the ambiguity in the solution can 
be understood as a rigid rotation of the Poincare sphere, 
which leaves the leakage solutions unchanged relative to 
each other. 
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Appendix A: Non-uniqueness of solutions 

We discuss the non-uniqueness of the solutions of Eq. 
and possible convenient conventions for choosing a specific 
solution. One obvious degeneracy is that multiplication 
of all the a's by one common phase factor independent 
of antenna, and all the g's, by another, in general differ- 
ent, common factor, does not affect the right hand side 
of Eq. |^. Also, the equation was written in a specific ba- 
sis, say right and left circular. But it would have had the 
same form when using any other orthogonal pair as basis 
since the source is unpolarized. Hence we are free to apply 
this change of basis to one solution to get another solution 
of Eq. ||. Under such a change, the coefficients transform 
according to 




(A.l) 

It is easy to verify that under this change, ot^a* 1 + 
g[g*' — ciicv* + gig]- Clearly, since x 2 1S unchanged by 
these transformations, an iterative algorithm will simply 
pick one member of the set of possible solutions, deter- 
mined by the initial conditions. Having found one such, 
one could apply a suitable transformation to obtain a so- 
lution satisfying some desired condition. For example, if 
one has nominally linear feeds, one might impose the sta- 
tistical condition that there is some mean linear basis with 
respect to which the leakage coefficients will be as small as 
possible. Such a condition has the advantage that a perfect 
set of feeds is not described in a roundabout way as a set of 
leaky feeds with identical coefficients, simply because the 
basis chosen was different. Carrying out the minimization 



of \9i\ 2 by the method of Lagrange multipliers, subject 
to a constant x 2 , we obtain the condition that ^2 a *9i = 0- 
This solution can be interpreted as requiring the leakage 
coefficients to be orthogonal to the gains, and is reasonable 
when we think about the opposite kind of situation, when 
the leakages are "parallel" to the gains, i.e. identical apart 
from a multiplicative constant. In such a case, we would 
obviously change the basis to make the new leakage zero. 
If we have a solution which does not satisfy this orthogo- 
nality condition, we can bring it about in two steps. First, 
choose an overall phase for the a's so that ^ ct*gi is real. 
Then, carry out a rotation in the g — a plane by an angle 
satisfying tan 6 — a *9i/(J2{9i9* ~ c^a*). This rotation 
has been so chosen that it makes the leakage "orthogo- 
nal" to the gains, in the sense required above. Even after 
this is done, we still have the freedom to define the phase 
zero independently for the two orthogonal states. This is 
because we are only dealing with unpolarized sources. Of 
course, if we had a linearly polarized calibrator, the rela- 
tive phase of right and left circular signals would not be 
arbitrary. 

A more geometric view of this degeneracy is obtained 
when we think in terms of the Poincare sphere represen- 
tation of the states of polarization of all the feeds. The 
cross correlation between the outputs of two feeds, both 
of which receive unpolarized radiation, has a magnitude 
equal to the cosine of half the angle between the represen- 
tative points on the sphere. Measurements of all such cross 
correlations with unpolarized radiation fixes the relative 
geometry of the points on the sphere, while leaving a two 
parameter degeneracy corresponding to overall rigid rota- 
tions of the sphere. This degeneracy can be lifted by the 
measurement of one polarized source at many parallactic 
angles. 

Finally, we note that for the purpose of correcting the 
observations of unpolarized sources for the effects of non- 
identical feed polarization, the degeneracy is unimportant, 
because the correction factor is precisely the right hand 
side of Eq. || which is unaffected by all the transformations 
we have discussed. 
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